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I The relaxation of hot nuclear matter to an equilibrated state in the central zone of heavy-ion 

' collisions at energies from AGS to RHIC is studied within the microscopic UrQMD model. It is 

' found that the system reaches the (quasi) equilibrium stage for the period of 10-15 fm/c. Within 

I . this time the matter in the cell expands nearly isentropically with the entropy to baryon ratio 

O I S/A = 150 — 170. Thermodynamic characteristics of the system at AGS and at SPS energies at 

^ I the endpoints of this stage are very close to the parameters of chemical and thermal freeze-out 

' extracted from the thermal fit to experimental data. Predictions are made for the full RHIC 

' energy ^/s = 200 AGeV. The formation of a resonance-rich state at RHIC energies is discussed. 

rS ; 

. Thermalization and chemical equilibration of hot and dense nuclear matter produced in 
ultrarelativistic heavy-ion collisions is a topic of great importance for the interpretation 
of current SPS and RHIC results. In our investigation we study the relaxation process in 
central Au-|-Au collisions at ^/s = 200 AGeV within the microscopic transport UrQMD 
model 0. Earlier studies at AGS and SPS energies revealed that the central reaction 
volume defined by a cubic cell of volume V = 125 fm^ is well suited for this kind of study 
01 . It contains enough particles to be treated as a statistical system, and its macroscopic 
characteristics become isotropic after some time. Figure |I] depicts the velocity distribu- 
tions of hadrons in the cell in transverse {x and y) and in longitudinal (z) directions. 
At t = 3 fm/c the longitudinal velocity distribution differs considerably from the distri- 
butions in the transverse plane, while at t = 5 fm/c the magnitudes and widths of all 
three distributions become very close to each other. Isotropy of the velocity distributions 
results in the isotropy of pressure in the cell. Pressure in longitudinal and in transverse 
direction is shown in Fig. H(a) for AGS, SPS, and RHIC energies. It is widely believed 
that the thermalization at RHIC sets in quite early. Indeed, at RHIC the pressure in the 
cell becomes isotropic at t ^ 5 fm/c compared with t ^ 8 fm/c (SPS) and t ^ 10 fm/c 
(AGS). Starting from t = 5 fm/c the results of the microscopic calculations are compared 
with the predictions of the statistical model (SM) of an ideal hadron gas 0. The val- 
ues of the energy density e, baryon density pb, and strangeness density ps, determined 
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microscopically, are used as an input to obtain particle yields, partial energy densities, 
pressure, and entropy density via the temperature T, baryocliemical potential fiB, and 
strangeness chemical potential fis- As seen in Fig. |](a), the microscopic pressure is very 
close to the grand canonical pressure after the onset of the (quasi) equilibrium stage. 



300 
200 

N 

>> 

^^100 

II 

"l 
150 

■« 100 

50 





t= 3 fm/c 



4- 
4- -f 



o - v.. 



ft ♦ ♦ 



M | 

t= 5 fm/c 



:§ 10-' 




-0.5 



0.5 




Au+Au Vs=200 AGeV 



Pb+Pb 1 60 AGeV 



Au+Au 10.7 AGeV 



t (fm/c) 



Figure 1. Hadron velocity distributions 
dN/dvi, i = X (□), y (Q)), and z (•) at 
t = 3 fm/c and t = 5 fm/c in a central 
cell in Au+Au collisions at RHIC. 



Figure 2. Time evolution of the components of 
the pressure tensor P (a) and of the entropy per 
baryon S/A (b) in the central cell of heavy ion 
collisions compared to the SM results. 



The entropy density per baryon S/A, defined both microscopically and macroscopically 
(via the Gibbs thermodynamic identity Ts = e + P — ^bPb — l^sPs) is presented in 
Fig. ^(b). The expansion in the cell proceeds nearly isentropically with S/A = 12 (AGS), 
32 (SPS), and 150 (RHIC) (cf. (S/A)*'^^™-^* ^ 14 (AGS) and = 36 (SPS) |]). Therefore, 
it would be very interesting to compare the UrQMD estimate (S'/A)^^^*-' = 150 — 170 with 
the value extracted from the thermal model fit to RHIC experimental data. 

The equation of state (EOS) in the (P, £:)-plane is shown in Fig. ^ For all three 
energies it can be well approximated by a simple linear dependence P/e = 0.12 (AGS), 
and 0.15 (SPS and RHIC). Note, that this version of the model does not imply the 
formation of the quark-gluon plasma, therefore, there are no kinks in this plot that can 
be attributed to quark-hadron phase transition. The evolution of the EOS in the (T, /xs)- 
plane is depicted in Fig. | together with the chemical freeze-out and the thermal freeze-out 
parameters obtained from a thermal fit to experimental data at AGS and SPS energies. 
One can see that temperatures and chemical potentials in the cells at the beginning 
and at the end of the equilibrated stage are close to the thermodynamic parameters of 
chemical and thermal freeze-out, respectively. The UrQMD predicts that at chemical 

195 ± 5 MeV and /i™FO = 43 ± 2 MeV, while at thermal freeze- 
130 ± 5 MeV and /i^™m.FO = 43 ± 3 MeV, i.e., the evolution of nuclear 
matter in the cell proceeds at constant baryon chemical potential. Calculations show 
very weak (within few MeV) difference between the results for Au+Au at y/s = 130 
AGeV and at ^/s = 200 AGeV. A thermal fit to particle ratios measured by the STAR 
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Collaboration (preliminary data) yields Tchem.FO = 190 ±20 MeV and /ig = 15±5 MeV 0, 
i.e. hb = 45 ± 15 MeV, which is in remarkable agreement with the UrQMD calculations. 
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Figure 4. The same as Fig. 3 but for the (T, //b)- 
plane. Sohd symbols correspond to the stage of ki- 
netic equilibrium, open symbols indicate the pree- 
quilibrium stage. Symbols with error bars are chemi- 
cal and thermal freeze-out parameters obtained from 
the SM fit Q. The hatched area shows the expected 
region of the quark-hadron phase transition. 

It is interesting to check the correspondence of the cell conditions at chemical freeze-out 
to the criterion {E)/{Nh) ~ (Mh) + 3/2T ^ 1 GeV introduced in Q]. The mean energy 



Figure 3. The evolution of pressure P 
and baryon density £ in the central cells 
of the heavy- ion collisions at AGS, SPS, 
and RHIC energies. 



Table 1 

Energy per hadron {E)/{Nh), energy density 
e, and temperature T in the cell at the be- 
ginning of kinetic equilibrium. 





time 


e 


{E)/{N,) 


T 




fm/c 


GeV/fm3 


GeV 


MeV 


AGS 


10 


0.68 


1.11 


129 


AGS 


12 


0.39 


1.06 


116 


SPS 


8 


0.74 


0.88 


170 


SPS 


10 


0.46 


0.80 


161 


RHIC 


5 


2.33 


1.08 


201 


RHIC 


6 


1.70 


1.01 


193 



per particle at the beginning of the equi- 
libration in the cell at AGS, SPS, and 
RHIC energies is listed in Table |l[ With 
rising bombarding energy from AGS to 
SPS there is a transition from baryon to 
meson dominated matter. It leads to the 
drop of (M/j) from nucleon mass to mass 
of p meson in accord with With fur- 
ther increase of the freeze-out tempera- 
ture the yields of heavy meson resonances 
rise faster than that of light mesons, thus 
leading to the rise of (M/j) with y/s. 



This means that not only the temperature, but also the mean mass of a particle is in- 
creased in the cell at RHIC energies, i.e., there should be more heavy resonances compared 
with the SPS cell. Therefore, the ratios of hadronic abundances are studied (see Fig. |). 
Here the results are presented separately for non-strange and strange baryons and mesons. 
In the baryon sector the resonances dominate over the strange and non-strange baryons 
until the end of the simulations. This can be taken as an indication of the creation of 
long-lived resonance-rich matter. The fraction of baryon resonances is almost 70% of all 
baryons in the cell at RHIC at 5 < t < 19 fm/c, while at SPS and AGS the number of 
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baryon resonances decreases from 70% to 35%, and from 60% to 25%, respectively. The 
meson fractions of resonances shrink within the time interval 5 < t < 19 fm/c from 60% 
to 30% (RHIC), 50% to 20% (SPS), and 40% to 15% (AGS). But at RHIC energies the 

hot hadronic matter in the cell as well as in 
the whole volume of the reaction is meson 
dominated. The mesons, baryons, and an- 
tibaryons carry 90%, 7%, and 3% of the to- 
tal number of particles in the RHIC cell at 
t > 10 fm/c (cf. 85%, 14.5%, 0.5% at SPS 
and 50%, 50%, 0% at AGS). The microscopic 
ratios for mesons (Fig. ^, right panels) seem 
to be very close to the SM ratios. Since the 
freeze-out occurs at t ~ 21 fm/c in the cen- 
tral cell at RHIC energies, the matter in the 
cell is frozen before reaching complete chem- 
ical equilibrium. This circumstance compli- 
cates the extraction of the chemical and ther- 
mal freeze-out parameters by means of the 
standard thermal model fit. 
The rapidity distributions of baryon reso- 
nances in Au-|-Au collisions at ^/s = 200 
AGeV are found to be nearly flat in the 
rapidity interval \y\ < 3.5 0. More than 
80% of the baryon non-strange resonances 
are still A's (1232). The density of directly 
reconstructible baryon resonances, that de- 
cay into final state hadrons, per unit rapid- 
ity at RHIC is quite high, and the resonance 
rich matter can be detected experimentally. 
The results of our study can be summarised as follows. The formation of long lived 
resonance- abundant matter is found. UrQMD predicts that T = 195 ±5 MeV, fiB = 43 ±2 
MeV, 150 < S/A < 170, and {E)/{Nh) ^ 1 GeV at chemical freeze-out in central Au+Au 
collisions at RHIC. The equation of state has a linear dependence P = 0.15 e. The 
UrQMD cell calculations show that strangeness to entropy ratio monotonically increases 
with rising ^/^ as Rs = 0.025 (AGS), 0.04 (SPS), and 0.05 (RHIC). 
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Figure 5. Time evolution of the hadron 
to resonance ratio R = H/Hj-es in 
the central cell of Au-|-Au collisions at 
RHIC shown separately for baryons, an- 
tibaryons, and mesons as well as for non- 
strange hadrons, strange hadrons, and 
total hadron yields. Circles denote the 
UrQMD predictions, stars correspond to 
the SM results. 



